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Abstract 

In this paper we present a mathematical analysis of the photoelectric effect for one- 
electron atoms in the framework of non-relativistic QED. We treat photo-ionization as 
a scattering process where in the remote past an atom in its ground state is targeted 



|. by one or several photons, while in the distant future the atom is ionized and the elec- 



tron escapes to spacial infinity. Our main result shows that the ionization probability, 
to leading order in the fine-structure constant, a, is correctly given by formal time- 
dependent perturbation theory, and, moreover, that the dipole approximation produces 
an error of only sub-leading order in a. In this sense, the dipole approximation is 
rigorously justified. 

> , 

q ■ 1 Introduction 

00 

Even today, more than 100 years after its discovery by Hertz, Hallwachs and Lenard, the 
phenomenon of photoionization is still investigated, both experimentally and theoretically 
' [HII]- This research is driven by novel experimental techniques that allow for the production 

of very strong and ultrashort laser pulses. In contrast, the photo electric effect in the early 
experiments is produced by weak, non-coherent radiation of high frequency. There is a third 
^ . physical regime, where the radiation is weak, of high frequency, and coherent. This regime is 

the subject of the present paper. We consider one-electron atoms within the standard model 
of non-relativistic QED, and we present a mathematically rigorous analysis of the ionization 
process caused by the impact of finitely many photons. Improving on earlier results con- 
cerning more simplified models, we show that the probability of ionization, to leading order 
in the fine-structure constant, is proportional to the number of photons, and, in the case of 
a single photon, it is given correctly by the rules of formal (time-dependent) perturbation 
theory. It turns out that the dipole approximation produces an error of subleading order, 
which provides a rigorous justification of this popular approximation. 

Let's briefly recall the standard model of one-electron atoms within non-relativistic 
QED. More elaborate descriptions may be found elsewhere [20|, 13 lj . States of arbitrarily 
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many transversal photons are described by vectors in the symmetric Fock space 

^:=0<Sn[® n L 2 (M 3 x{l,2})] 

n>0 

over L 2 (M 3 x {1, 2}). Here S n denotes the projection of L 2 (M 3 x {1, 2}) n onto the subspace 
of all symmetric functions of (ki, Ai), . . . , (k n , A n ) G M 3 x {1, 2}, and S Q L 2 (M. 3 x {1, 2}) := C. 
We shall use Q, to denote the vacuum vector (1,0,...) G T. Nf is the number operator 
in J 7 , and Hf = dT(uj) denotes the second quantization of multiplication with uj{k) = \k\ 
in L 2 (IR 3 x {1,2}). See [27] . X.7, for the notation dr(-) and for an introduction to second 
quantization. The creation and annihilation operators a*(h) and a(h), for h G L 2 (M 3 x 
{1,2}), are densely defined, closed operators with a*(h) = a(h)* and with 

[a*(/t)*]W = ^S n {h ® *("- 1 )) 

for vectors ^ = , . . .) from the subspace D(Nj ). Here, \lK n ) denotes the n- 

photon component of 

The system studied in this paper is composed of a non-relativistic, (spinless) quantum 
mechanical, charged particle (the electron), and the quantized radiation field which is cou- 
pled to the electron by minimal substitution. In addition, there is an external potential V, 
which may be due to a static nucleus. The Hilbert space is thus the tensor product 

H := L 2 {R 3 ) 0.F, 

and the Hamiltonian is of the form 

H a = (p + a§A(ax)) 2 + V + H f (1.1) 
= H + W, 

where Hq = H e \ + Hf, H e \ = — A + V, and W = H a — Ho. The quantized vector potential 
A(ax), for each x G M 3 , is a triple of self-adjoint operators, each of which is a sum of a 
creation and an annihilation operator. Explicitly, 

A(ax) = o(G x ) + a*(G x ), G x (k, A) := -^=s(k, A)e" iak - X , (1.2) 

V 2|k| 

where e(k, A) G M 3 , A = 1,2, are orthonormal polarization vectors perpendicular to k, and 
k is an ultraviolet cutoff chosen from the space 5(M 3 ) of rapidly decreasing functions. No 
infrared cutoff is needed. Here and henceforth, the position of the electron, x G M 3 , and 
the wave vector of a photon, k 6 I 3 , are dimensionless and related to the corresponding 
dimensionfull quantities X, K by X = (ao/2)x and K = (2a/ao)k, where ao := h 2 /me 2 
is the Bohr-radius, m > is the mass of the particle, e its charge, and a = e 2 /he is the 
fine structure constant. It follows that X • K = ax • k, and in units where h, c, and four 
times the Rydberg energy 2ma 2 are equal to unity, the Hamiltonian of a one-electron atom 
with static nucleus at the origin takes the form with V(x) = —Z/\x\, Z being the 

atomic number of the nucleus. For simplicity, we confine ourselves, in this introduction, to 
this particular potential V. In nature, a ~ 1/137, but in this paper a is treated as a free 
parameter that can assume any non-negative value. 
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For all a > 0, the Hamiltonian H a is self-adjoint on D(Hq) and its spectrum o~(H a ) is 
a half-axis [E a ,oo) [Ml [23]. Moreover, 

E a := inf a(H a ) 

is an eigenvalue of H a , and, at least for a sufficiently small, this eigenvalue is simple [2l [2Tj . 
We use <& a to denote a normalized eigenvector associated with E a . Another important 
point in the spectrum of H a is the ionization threshold S a , which, for our system, is given 
by S Q = inf o~(H a — V). In a state vector from the spectral subspace Ranl^^ ^) 
the electron is exponentially localized in the sense that 

||e^l ai ll ( _ 00>Ea _ e] (H a )||<oo (1.3) 

for all with j3 2 < e PS|. 

The phenomenon of photo-ionization can be considered as a scattering process, where 
in the limit t — > — oo, the atom in its ground state is targeted by a (finite) number of 
asymptotically free photons, while in the limit t — > oo the atom is ionized in a sense to be 
made precise. We begin by discussing incoming scattering states and their properties. To 
this end it is convenient to introduce the space L^j(M. x {1, 2}) of all those / G L 2 (R x {1,2}) 
for which 

11/11*:= £ / |/(k,A)| 3 (l + a;(k)- 1 )rf , A ! <oo. (1.4) 

A=l,2 J 

Given / G L^(R x {1,2}), the asymptotic creation operator a* (/) is defined by 

a* (/)*:= lim e iH ° t a*(f t )e- iH " t V, ft := e'^f, (1.5) 

t— oo 

and its domain is the space of all vectors Vl/ G Dd^ffol 1 / 2 ) for which the limit (jl.5j) exists. 
This is know to be the case, e.g., for the ground state \& = f» a . Moreover, it is known that 
a -(/i) • • • a*(/n)^a is well defined and that 

e-^aK/i) •••al(/ n )# Q 
= a*(/ 1 , i )---a*(/ n , t )e- i ^^ Q + (l), (t - -oo), (1.6) 

whenever G ^(lx {1,2}) for alH = 1, . . . , n [22]. By (JTSJ, a* (/ a ) • • • a*_(f n )$ a 

describes a scattering state, which, in the limit t — > — oo is composed of the atom in its 
ground state and n asymptotically free photons with wave functions fi,...,f n . Results 
analogous to those on o!!L(/) hold true for the asymptotic annihilation operators a~(/) |22j . 

The asymptotic annihilation and creation operators satisfy the usual canonical commu- 
tation relations: e.g. 

[a-(f),at(g)] = (f,g) (1.7) 

for all f,g G L 2 (M 3 x {1,2}). Moreover, the ground state $ Q is a vacuum vector for 
asymptotic annihilation operators in the sense that 

a_(/)$ a = for all / G L 2 (M 3 ). (1.8) 
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Hence, if / = (fi,m 1 , . . . , f n , m n ) £ [L 2 (R 3 x {1, 2}) xN]" with (f h fj) = (%, then it follows 
from ([LTD and ([L8]) that 

n 

is a normalized vector in "H. All these properties of a_(/), a* (/) hold mutatis mutandis 
for the asymptotic operators a+(g), defined in terms of the limit t — > +oo. 

We are interested in the probability that e~ lHat a*_(f)<& a describes an ionized atom in 
the distant future, but we are not interested in the asymptotic state of the electron or 
the radiation field in the limit t — > +oo. We therefore shall not attempt to construct 
outgoing scattering states describing an ionized atom, which is a difficult open problem. 
Instead we base our computation of the probability of ionization on the following reasonable 
assumption: the atom described by e~ lHat a*_(f)<& a is either ionized in the limit t — > do, or 
else, in that limit, it relaxes to the ground state in the sense that e~ iHat a*_{f)^ a , for t large 
enough, is well approximated by a linear combination of vectors of the form 

a*(gi,t)...a*(g n>t )e- iE ^ a . (1.10) 

More precisely, relaxation to the ground state occurs if a*L(f)$> a belongs to the closure of 
the span of all vectors of the form 

a* + (gi) ■ ■ ■ al(g n )<!> a = lim e iHat a*(g x>t ) . . . a* (g n>t )e- iEat ^ a , 

with gi,ujgi £ L^(IR 3 x {1,2}). Let H+ denote this space and let P? be the orthogonal 
projection onto Then ||P°al(/)$ Q || 2 is the probability for relaxation to the ground 
state and 

1 - ||P>1(/)<D Q || 2 = ||(1 - P>l(/)<D a || 2 (1.11) 

is the probability of ionization. 

The assumption that relaxation to the ground state is the only alternative to ioniza- 
tion, is motivated by the conjecture of asymptotic completeness for Rayleigh scattering, 
which is the property, that every vector ^ £ Ti describing a bound state in the sense 
that sup t ||e £ ' x 'e _ *" Ha ' , I'|| < oo for some e > 0, will relax the ground state in the limit 
t — > oo. In view of (II. 3p . asymptotic completeness for Rayleigh scattering implies that 
1(— oo,s a ) C^a)) which can be proven for simplified models of atoms [301 \E[ \TE[ I14j. 

The following two theorems will allow us to compute (jl.lip . 

Theorem 1.1. Suppose that f±, . . . , f n £ L 2 (R 3 x {1, 2}) where Yl1=i ^)/»(") ^) belongs 
to Cg(R 3 \{0}, C 3 ) for each i, and let / = (ft, . . . , f n ). Then: 

/oo 
2p(s) ¥ > ei ® [A(0,s),a*(l)]nds + O(a 5 / 2 ) (1.12) 
-oo 

where p(s) = e iH -' s pe~ iH - lS and A(0,s) = e iH f s A(0)e- iH f s . 
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The first term of (jl. 12[) gives no contribution to the ionization probability (jl.lip because 
a+(f)$> a £ H%- The second term is proportional to a 3 / 2 and it is due to scattering processes 
where one of the n photons /i, . . . , f n is absorbed. The remainder terms are of order 0(a 5 / 2 ) 
and stem from the dipole approximation A(ax) — > A(0), from dropping a 3 A(x) 2 and from 
ignoring processes of higher order in a 3 / 2 . To isolate the contribution of order o? from 
(fTTTT]) using CCS]), we need: 

Theorem 1.2. Suppose that7i.+ D l(-oo,s a )(^a) for a in a neighborhood ofO, and suppose 
that H e \ has only negative eigenvalues. Then 

lim lP? = lppCfld) (1.13) 

a— >U 

in £/ie strong operator topology. 

Combining Theorem 1 1 . 1 1 and Theorem 11.21 we see that 

||(1 - P«)a*_(/)$ Q || 2 = ||(1 - F?)(a* (/) - a;(/))$ a || 2 

= ||l c (i^ e i) («K/) " <(/)) $ a || 2 + o(a 3 ) 

/oo 
p(«) VeJ 8) [A(0,s),a*(/)]ttds|| 2 + o(a 3 ) 
-oo 

where l c (iT e i) = 1— l pp (i/ e i), and where the second equation is justified by the a dependence 
of a*_{f)<& a — a+(f)$ a as given by (|1.12|) . We are now going to express the coefficient of 
a 3 in terms of generalized eigenfunctions of H c \, which makes it explicitly computable in 
simple cases. A general and sufficient condition for the existence of a complete set of 
generalized eigenfunctions is the existence and completeness of a (modified) wave operator 
0+ associated with H e \. This condition is satisfied for our choice of V. It means that there 
exists an isometric operator Q + G C(Tt e i) with Ran$7 + = l c (H e \)TC e i and H C \Q + = $7 + (— A). 
In particular, the singular continuous spectrum of H e \ is empty. Given the wave operator 
fi+ and the fact that (H c \ — i)~ 1 (x)~ 2 is a Hilbert-Schmidt operator, it is easy to establish 
existence of generalized eigenfunctions ^> q , q G M 3 , of H e \ with the following properties |26j : 

(i) The function (x, q) (q)~ 2 (x)~ 2 <y9 q (x) is square integrable on R 3 x M 3 , in particular 
(x)-V q G L 2 (M 3 ) for almost every q G R 3 . (x) := (1 + (x) 2 ) 1 / 2 . 

(ii) If V G £>(|x| 2 ) then 

lllacGEWII 2 = / \(^^)\ 2 d 3 q (1.14) 

(iii) If F : R -» C is a Borel function, V G Z^(|x| 2 ) n D(F(H c i)), and F(H c \)ip G £>(|x| 2 ), 
then 

(<p v F(HM =F(q 2 )<<^> (1-15) 

for almost every q£l 3 . 

In (ii) and (iii) we use (ipq,ip) to denote the integral J y> q (x)^(x) d 3 x, which is well defined 
by (i) and by the assumption ip G L>(|x| 2 ). 

The Theorem 1 1.1 1 in conjunction with (i)-(iii) implies the following theorem, which is our 
main result specialized to the case of only one asymptotic photon in the incident scattering 
state. 
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Theorem 1.3. For all f £ L 2 (R 3 x {1,2}) with Ea=i A )/(-> A ) G C^(M 3 \{0}, C 3 ), 



\l ac (H el )(at(f)$ a -a* + (f)$ c 



a 3 / d 



((Pq, ynp e l) 



e i(q 2 -E )t (n,E(0,t)o*(/)fi) dt 



+ 0{a A ) (1.16) 



cis q — ► 0. Here, E(0,i) = — A(0, t)], <p e i is a normalized ground state of H c \ and <pn, 
q 6 M 3 , is any family of generalized eigenfunction of H c \ with properties (i)-(iii) above. 



The expression (|1.16p for the ionization probability can be understood, on a formal level, 
by first order, time-dependent perturbation theory. To this end one considers the transitions 
'-Pel <S> / > (/? q (8> 0, for fixed q G R 3 , in the interaction picture defined by Hq. Then the 
time-evolution of state vectors is generated by the time-dependent interaction operator 
W{t) = e iHot We- iHot = 2a 3 / 2 p(t) • A(ax,t) + a 3 A(ax,t) 2 with p(i) = e iH e it pe -iH el t and 
A(ax, t) = e lHot A(a~x.)e~ lHot . In the computation of the transition amplitude to the order 
o?l 2 one drops a 3 A(ax, t) 2 and one replaces A(ax, t) by A(0, t), which is known as the 
dipole approximation. Then, an integration by parts using that 



2p(t) = |x(t), 



^A(0,t) =E(0,t), 



leads to a result for the transition amplitude which agrees with the expression in (|1.16|) 
whose modulus squared is integrated over q£K 3 . The Theorem 11.31 and its proof justify 
this formal derivation and the use of the dipole approximation. Note that ax = X, hence 
the ionization probability is of order a 3 rather than of order a, as a formal computation, 
similar to the one above, in dimension-full quantities would suggest. 

We prove a more general result than Theorem 11.31 where the incoming scattering state 
may contain several asymptotic photons, and where the external potential V is taken from a 
large class of long range potentials. In the case where the asymptotic state at t = — oo is of 
the form (jl.9p and each of the photons fi, . . . , f n G L 2 (M 3 x {1, 2}) satisfies the hypotheses 
of Theorem 11.31 in addition to (fi, fj) = 5{j, our result says that 

n 

||l ac ( J ffei)(aI(/)$a-<(/)*a)|| 2 = a 3 ^m z pW(/ z ) + C»(a 4 ) (1.17) 



l=i 



with 



P {3 \fi) 



d 3 q 



/oo 
<n,E(o,t)o*(/,)n> e'tf-^dt 
-oo 



(1.18) 



The integral with respect to t in (|1.18p can be computed explicitly in terms of fi and Go, 
and it gives 



e i{^-E )t (o,E(0,i)a*(/i)fi) dt 

K(k)^2\k\ £ e A (k,A)/ l (k,A)d<r(k) 1 



ITT 



k|=q2-E 



A=l,2 



G 



where dcr(k) is the surface measure of the sphere {k G R 3 : |k| = q 2 - E } in R 3 . The 
integration over the spheres with |k| = q 2 — Eq expresses the conservation of energy in the 
scattering process, and the additivity (|1.17p of the ionization probability with respect to the 
incoming photons corresponds to the experimental fact, that the number of photo-electrons 
is proportional to the intensity of the incoming radiation. 

In Section Owe give a second derivation of a 3 Ps(f) based on a space-time analysis of 
the ionization process. This approach, in a slightly different form, was introduced in the 
papers [31 [33], and does not assume asymptotic completeness of Rayleigh scattering. 

The existence of outgoing scattering state describing an ionized atom and an electron 
escaping to spacial infinity is a difficult open problem in the model described above. Only 
for V = such states have been constructed so far [251 E]- Hence it is not possible yet 
to study the ionization probability based on transition probabilities between asymptotic 
states. 

Previously ionization by quantized fields was investigated in [31 HH [T71 [33] . [3] and [33] 
are precursors of the present paper on simpler models of atoms and the ionization prob- 
ability defined in a different, but equivalent way. In [16|, [T7] it is shown that a thermal 
quantized field leads to ionization in the sense of absence of an equilibrium state of atom 
and field. 

There is a large host of mathematical results on ionization by classical electric fields: 
Schrader and various coauthors study the phenomenon of stabilization by providing upper 
and lower bounds on the ionization probability, see [10|, IT2l [TT] and the references therein. 
They use the Stark-Hamiltonian with a time dependent electric field £{t) that vanishes 
unless < t < r < oo. Lebowitz and various coauthors compute the probability of ioniza- 
tion by an electric field that is periodic in time; see [29] and references therein. Most of 
these papers study one-dimensional Schrodinger operators with a single bound state that 
is produced by a <5-potential. Ionization in a three-dimensional model with a <5-potential is 
studied in [6]. 

Acknowledgment: M.G. thanks Vadim Kostrykin for pointing out that it is advantageous 
to define ionization as the opposite of binding. 

2 Notations and Hypotheses 

For easy reference, we collect in this section the definitions, our notations and all hypotheses. 
As usual, L 2 (M 3 x {1, 2}) denotes the space of square integrable functions / : K 3 x {1, 2} — > C 
with inner product 



We recall from the introduction that L 2 (M 3 x {1, 2}) consists of those functions / E L 2 (M 3 x 
{1,2}) for which the norm H/H^ defined in (II ,4|) is finite. Regularity assumptions will be 
imposed on the vector-valued function 
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(2.1) 



A=l 
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rather than on on /(•, 1) and /(•, 2). It is useless to impose smoothness conditions on /(•, A) 
because it is (|2.ip that matters and because the polarization vectors e(k, 1) and e(k, 2) are 
necessarily discontinuous. On the other hand, every square integrable function / : M 3 — » C 3 
with k • /(k), for a.e. k € M 3 , can be approximated, in the L 2 -sense, by smooth functions 
of the form (|ZTJ) . 

It is convenient to collect a family /i, /jv G L 2 (M 3 x {1,2}) of photon wave functions 
in an 7V-tupel / = (/i, /jv)- We define 

a(/) := a^O-.-a^) 
a*(/) := a*(/i) •■•£**(/*). 

This should not lead to confusion with (|1.9jl . where / also includes occupation numbers. 
For the various parts of the interaction operator W = H a — Hq, we use the notations 

W dip := 2p-A(0), 
W {1) := 2p-A(ax), 
W {2) := A(ax) 2 . 

It follows that 

W = Jw^ + a 3 W {2) = J\V dip + 0(ai ) 

where the last equation is purely formal, but we shall give it a rigorous meaning in this 
paper. The Hamiltonian 

H a = H + W 

is self-adjoint on the domain of — A+Hf provided that V is infinitesimally operator bounded 
with respect to —A, [23, [23]. This is the case, e.g., if V is the sum of Coulomb potentials 
due to static nuclei; all our results are valid for such V. Nonetheless, it is useful to identify 
the properties of V that are essential for our analysis. From now on, we shall only assume 
the following hypotheses on V: 

Hypotheses: Both V and ~x.-W belong to G L 2 oc (IR 3 ), \\m^^ oc V{ii) = 0, and there exist 
constants fi > and R > such that for \j3\ = 1,2 we have 

\d?V(x)\ < Ixr^l-", if\x\>R. 

Moreover, E := inf a(H e \) < 0. We define e\ := mi(a(H c \)\{E }). 

From these Hypotheses it follows that a ess (H e \) = [0, oo), that a sc (H c \) = and that Eq 
is a simple eigenvalue. In fact, the decay assumptions on V imply long-range asymptotic 
completeness |9J , which is what we use to infer the existence of a complete set of generalized 
eigenfunctions. All this remains true if a singular short-range potential is added to H e \. 

The time evolution of an operator B in the interaction picture will be denoted by B(t), 
that is, 

B(t) := e iHot Be- iHot , 

and B t := B(-t). Note that p(t) = e iH ^pe- iH ^ , A(0,t) = e iH f l A(0)e~ iH f l and that 
a#H t ) = e- iHot a*(f)e iHot = a#(f) t . 
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3 Commutator estimates and scattering states 

The main purpose of this section is to establish bounds on the commutators [W^\a*(f )] 
applied to $ a for G {W^\ W^ 2 \ W^}. We are interested in the decay as |t| — ► oo 

and in the dependence on a. Typically, our estimates are valid for a < a, where a is 
defined in Proposition IA.31 As a simple application of our decay estimates in t, we will 
obtain existence of the scattering states 

4(/)$ Q = lim e ltH «a*(l t )e- UH «<t> a , 

which was already established in [22] in larger generality. Here f t = (/i,t 3 /jv,i) an d 
f jt ■= e- itw fj. Given I G {1, N}, we write 

«*(/[*],*) : = a*(/i,t)---a*(/ I -i, i )A(ax)o*(/ I+1 , t )---a*(/j V , t ) J 
a *(/(/),t) := **Ui,t)---a*Ui-i,t)a*Ul+i,t)---a*{fN,t). 



For x G K 3 , (x) := (1 + |x 



2 U/2 



Lemma 3.1. Suppose that f G L 2 (M 3 x {1,2}) e/ G C^(R 3 \{0}, C 3 ) for a given 

n G N. T/ien £/iere exists is a constant c\^ n = ci jn (/) smc/i i/iai 



I (G ,/ t ) | < ci, n 
|(G XJ /t)| < ci, n 
(G x -G ,/ t )| < ci, n 



1 



1 + 

l + (a|a;|) n 
1 + |t| n 

a|x|(ax) n 
1 + \t\ n 



for allxe K 3 , 
for all x G K 3 . 



(3.1) 
(3.2) 
(3.3) 



Proof. Estimate (]3.ip follows from ()3.2p . We next prove (13.20 . By a stationary phase 
analysis of 

K (k) 



(G x , ft, 



drk ' ; e- 

TMk) 



(e/)(k) 



(3.4) 



we obtain |(G x ,/ t )| < C n |t|- n for a|x| < |i|/2, [28] Theorem XI.14. It follows that 

C n 



(G X3 /t)|l{2a|x[<|t[} < j^p 
(G x , /t)|l{2a|x[>|i[} < C 



2a|x| 



where C := sup tgK xgK 3 |(G x ,/j)| < oo. This proves (|3.2|) , To prove (|3.3p we write 

(G x -G ,/ t ) = J e- u ^F x (k)d 3 k 



where 



F x (k) = iak ■ x 



«00 



(e/)(k) g (ak-x) 
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and g : R — ► C denotes the real-analytic function given by g(s) = (e ls — l)/(is) for s / 0. 
g and all its derivatives are bounded, and by assumption on /, F x G Co°(M 3 \{0}, C 3 ) for 
each x. It follows that 



sup 

x,keK 3 ,x^O 



<^F x (k) Ixl-^ax)-^! < oo, 



which implies (13.31) . again by stationary phase arguments. 



□ 



Lemma 3.2. Suppose that e/i,...,e/jv € Cg (R 3 \{0}, C 3 ) /or a given n G N, and let a be 
defined by Proposition \A.3[ Then there exist constants a > and C2 t n = C2,n(f), such that 
for alla<a,t£ R, and G {W^, W( 2 ), W dip }, 



(3.5) 



Proof. By definition of a*(f 



A' 



[W^\a*{[_ t )]<5> a = £ a*(/i,t) • • • «*(//-!,*) [^°'\ o*(/i,t)] «*(/m,t) ■ ■ ■ a*(/jv,*)$<* (3-6) 



i=i 



and by definition of and V^ 2 ) 

[W (1) ,a*C/U>] = 2(G x ,/ M )-p 
[W^\a*(f l>t )] = 2(G X ,/ M )-A(ax) 

From ([321), iCT r d32D, (USD and Lemma EH it follows that 

\\{WW,a*(fJ]$ a \\ < I ^-[|(ir / + l)^(ax>»p* 



+ . 

[W( 2 ),a*(/ 4 )]d> a || < _^-||(ff / + l)f( axr $ 



(3.7) 
(3.8) 



(3.9) 
(3.10) 



with some constant c n . Thanks to Lemma lA.4} these upper bounds are bounded uniformly 
in a < a, a being defined by Proposition IA.31 This proves (j3.5l) for j = 1, 2. The assertion 
for W di v now follows from W di P = W^\ x=0 , which leads to a bound for || [W di P, a*(/ t )]$ a || 
of the form (1331) with x = 0. □ 



Proposition 3.3. For all efi, e/jv G Cq (R 3 \{0}, C 3 ) i/iere exists a constant C3 = 03(f), 
such that for all a < a and for all s G R, 



W(i) -iy di P, a *(/ ) 



< 



c 3 a 



1 + s 2 ' 

Proof. By (13.60 for j = 1, (13. 7p . and the corresponding equations for 14^ dip 

AT 

_ ^ dip ,a*(/ s )]$ Q = 2^<G X - Go, /,,.) • pa* (I s )$ a 



(3.11) 



l=l 
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where 



(G x -G ,/i, a >-pa*(^ )* tt || < a^^|||x|(ax) 2 a*(/ m >* c 



1 + s 2 



< a^-^|||x| 2 (ax) 4 p$ a || 1/2 ||a(/ a)!s )a*(/ (0!e )p$ a || 1/2 



1 + s 2 



by (|3.3p and the Cauchy-Schwarz inequality. The norms in the last expression are bounded 
uniformly in a < a by Lemma I A. II and Lemma IA.4I □ 

Lemma 3.4. For allefi, ...,e/jv € Cq (R 3 \{0}, C 3 ), there exists a constant C4 = 04(f) < 00, 
such that for all a < a and s, t G R 



||[^P,a*(/ f )]($ a -$ )|| < 1+ C ,^ g|2 

||[< ip ,«*(Z,)]<i>o|| < 



1 + 1* -fi 

3 

04 a 



ll^^a^)]]*,,!! < 



(3.12) 
(3.13) 

(3.14) 



Proof. Since [W^ 113 , a*(fij)] = 2(Go, fij-s)-p.s, which commutes with the creation operators 

a*(fi,t), 



N 



[W^,a*(l t )] = Y, a *(f^)---a*(fi-i,t)[Wf^,a*(f l , t )]a*(f l 



l+i,t) 



a* a 



N,t) 



1=1 

N 



2 E a *(/mtK G o>/M- s >-P s 



(3.15) 



1=1 



where |(Go, fi,t-s)\ < q(1 + it — s ) 2 ) 1 by ()3.ip . In view of Lemma [A. II and Lemma [A. 51 
this proves (|3.12p . The proof of (13.130 is similar. 
From (I3.15P we obtain, that 



v 



W, 



1=1 



Hence, by (|3.1|) . it suffices to show that \\W ( - j ^a*(f^ t )p s $ a \\ and || t )p a W^* a || are 

bounded uniformly in t,s and a < a. We shall do this for a*(f, . )p s W^^ a only, the 
proofs in the other cases being similar. Let m > (N — l)/2. Then 

IK(/(o,*)p* w(1)< M 
3 

^ E IKC/W^/ + ir"Wi(#ci + o _1 (^ei + o(s> + ir^(ax)<D a || 
3=1 

3 

< C^||(F el + i)(F / + ir^(ax)«I> a || 
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with a constant C, that is finite by Lemma I A. 11 We now want to compare \\(H e \ + i)(Hf + 
l) m Aj(ax)$ a \\ with \\Aj(ax)(H el +i)(H f + l) m $ q,||, because the latter norm is bounded 
uniformly in a < a, by Lemma I A. II and by (|A,17p , Thus we compute the commutator of 
(H cl + i)(H f + l) m and A; (ax) = a*(G xJ ) + a(G xJ ) applied to $ a . Using 



3 



[H eh a*(G x j)] = a 2 a*(uj 2 Gy_j) - a } y 2a* (k m G X)j )p m 

m=l 

[(H f + ir,a*(G Xjj )] = J2{l) a *^ lG *MHf + l) m - 1 

and similar commutator equations for a(G K j), we see that all resulting terms have norms 
that are bounded, uniformly in a < a, thanks to (|A.17h and Lemma lA.ll □ 

For completeness of this paper we now use Lemma 13.21 to prove existence of the asymp- 
totic creation and annihilation operators on & a . More general results can be found in 

us E3- 



Proposition 3.5. Suppose / = (/i, ...,f N ) S [Ll(M 3 x {1,2})]^. Then, for all a < a, 

a* ± (f)* a := lim e^-VC/Je^-^a (3.16) 

exists, and 

||o±(/)*a|| <c 5 ||/iL---||/ivL, (3.17) 

wrai/i a constant c§ that is independent of a and f . Ifsfi £ Cq +1 (IR 3 \{0}, C 3 ) /orZ = 1, ...,N, 
then there exists a constant c n (f), such that 

c (f) 

||4(/)$ Q - e^a^Je-^*^!! < « 3/2 ^p (3-18) 
Proof. Suppose first, that efi,...,ef N G Co +1 (M 3 \{0}, C 3 ). Then 

^fe^a*(f )e-^$ Q ) = ie^-^pV, a* (/.)]*«, 



eft V 

and, by Lemma I3T21 



± 



[W^U^4d S <a^^ 



This estimate first proves existence of a±(f), by Cook's argument, and then it implies 
(ETTKl) . The existence of a±(/)$ a in the case where fj £ L 2 W (M 3 x {1, 2}) now follows from 
the approximation argument given in [22], Proposition 2.1. By the Lemmas IA.1I and I A. 41 

||e^«a*(/ t )e-^$ a || < \\a*{Q{H f + l)-T\\\\(H f + l)£$ a || 

< c 5 ||/x[|oj • • • II/aHL 

uniformly in i £ K and a £ [0, a]. Letting i — > ±oo in this estimate, we obtain (|3.17p . □ 
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4 Proofs of the main theorems 



4.1 A reduction formula 

In this section we first prove Theorem 14 . 1 1 b elow . which is a generalization of Theorem 1 
the latter corresponding to the choice r = 0. The generalization to arbitrary r £ R will be 
needed in Section El 

Theorem 4.1. lei efi, ef N E C£(R 3 \{0}, C 3 ). T/ten 
a;(/ r )$ a -a*(/ r )$ a 



3 



e-^ H °- Eo >2p(s)(p e i ® [A(0, s), a* (/)]« (is + ft(r, a) 



w/iere ||ft(r,a)|| = C(a 5 / 2 ) + C(a 3 |r|) as a -> 0. 

Remark. Part of the error C(a 5y/2 ) stems from passing to the dipole-approximation — > 
W dip . Hence its order 5/2 = 3/2 + 1 cannot be improved. 

Proof. Recall that Bt = B(—t) = e~ °Be °. To compare the time-evolutions generated 
by H a and Hq we will use that 

e i(H a -E a )t Bt $ Q = Bq?a + I e K Ha ~ E '* i '[%W,B a ]^ a da. (4.1) 

Jo 

This equation may be iterated because [iW, B s ] = [iW- s , B] s . From 

oi(/)* a = lim e^ H "- E «K*(f.)^ a 

— t— »±oo — 1 

and (|4.ip it follows that 

/oo 
-oo 

Only terms contributing to this integral of order c?l 2 need to be kept. Since 

a we may drop T-W 2 ), — W dip and restrict the interval of integration to |s| < a 1 

by Lemma 13.21 and Proposition 13.31 We obtain 

<(/)*«- a* (/)* Q 



io?' 2 



/oo 
e i(H a -E a )s [ W di P) a *(/J$ Q ds + 0(a 5/2 ) 
-oo 

f e i{Ha - Ea)s [W*Y,a*{r)} s <S> Q ds + 0(a 5/2 ). (4.2) 

./|s|<a-l 

Applying now (|4.1|) to the integrand in (|4.2|) and the time interval [r, s], rather than [0,s], 
we find 

I e< H -- E ->[W^,a*(f)} s 9 a ds 
= e* H "- E °> f [W^,a*(f)] T <S> a ds (4.3) 
+ [ ds f S H '- E ->[%W, [W?Z, a*(/ p )]]*„ dr. 
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By (13. 141) in Lemma 13.4} the norm of the double integral is bounded by 



const 



M + s 



|<a-i 1 + \ s \ 



q 3/2 ds = 0(a 3/2 \T\) + 0(a 3/2 ln(a)). 



(4.4) 



In the integral (|4.3f) we use Lemma 13.41 to replace $ a by $o an d to extend the integration 
over all set We find that 



\s\<a~ 1 



[W^,a*(f)] T <S> a ds 

[W^ s ,a*U T )]<S> ds + O{a) 

/oo 
e -t(fl- -B B )r[ W rdiP ) o *(/)]$ ^ + 0(a). 
-co 



oo 
oo 



(4.5) 



Equations (|43j) . P~5]h (fO|) and (|4"3j) prove the theorem because e" i( ^ Q "- E ° !) ' r a|(/)$ c 
o±(/ and because $o = Pel <8> ^- 

Theorem 14. II in the case r = becomes Theorem I l.li which implies that 
||l„ c (^ e i)(al(/)$ Q - <(/)$ a )|| 2 = a 3 pW(l) + 0(a 4 ) 

where 



□ 



l ac (F el ) / 2p(a) VcI ® [A(0,s),a*(/)]Ods 



P( 3 )(/) := 



We next show that P^ 3 \f) is additive in its one-photon contributions. 



(4.6) 



Proposition 4.2. Suppose that f = (/i,mi, . . . , f n , m n ) G [L 2 (M 3 ) xN]™ «;zf/i (/j, = 5^ 
and e/j G C 2 (M 3 \{0}, C 3 ). Then P$\f) = £f =1 ™iP {3) {fi) with 



P {3) (fi) 



oo 

l ac (H el ) J 2p(s)<p el -(n,A(0,s)a*(f l )n) ds 

— oo 
oo 

l ac (H cl ) J K(s)i Pel -{n,E{0,s)a*(f l )n) ds 



(4.7) 



Proof. Since a*(f) is a product of creation operators a*(fi) and since [A(0, s), «*(//)] 
(O, A(0, s)a* (fi)£l) , a scalar multiple of the identity operator, we have 

n 

[A(0, a), a*(/)] = ^ ^ <n, A(0, s)a*(/,)fi) a*(/,Jfi 



where = (fx, mi, ■ ■ ■ , fi, (mi — 1), ... , / n , m n ). The vectors a*(f^)Q are orthonormal 
by construction. Hence by definition of P^ 3 '(f) and by the Pythagoras identity, 

n 

p(3) (Z) = ^p(3) (//) 
1=1 
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with p( 3 \fi) given by the first equation in the statement of the proposition. The second 
equation in the proposition follows from 



2p(s)(p ei = —x(s)(p d , 



ds 



(O, A(0, s)a* (/,)fi) = - (O, E(0, s)a* (/,)«) 



by an integration by parts. The differentiability of s h- * x(s)ip e i and the expression for its 
derivative are established in Lemma |A. 61 □ 



4.2 Expansion in generalized eigenfunctions 

In this section we prove Theorem 11.31 and the stronger statement expressed by the Equa- 
tions (|1.17p and (|1.18|) . The ingredients are Theorem 11.11 Proposition 14.21 and a set of 
generalized eigenfunctions ipq with the properties (i)-(iii) in the introduction. Concerning 
the existence of y> q , we recall from [S], Theorem 4.7.1, that our hypotheses on V imply ex- 
istence and completeness of a (modified) wave operator fi + associated with H e \. Moreover, 
(H e \ — i)~ 1 (x)~ 2 is a Hilbert-Schmidt operator. 

Lemma 4.3. Suppose that <p : R — ► TL e i n D(|x| 2 ) is such that s i— ► <p(s) and s i— ► \x\ 2 ip(s) 
are continuous and absolutely integrable with respect to the norm ofH e i- Then ip(s)ds G 
£>(|x| 2 ) and 

r°° 2 

lac(H e \) / if(s)ds 



(<Pq,<p(s)) ds 



d A q. 



Proof. From the existence of the improper Riemann integrals ip(s)ds and \x\ 2 <p(s)ds 
and the fact that multiplication with |x| 2 is a closed operator, it follows that J^_ s ip(s)ds G 
D(\x\ 2 ) and that 



x 



(p(s)ds 



\x\ 2 ip(s)ds. 



This equation and property (i) of (fq imply that 



(p(s)ds 



W </>q ; 



\x\ 2 tp(s)ds 



/oo roc 
<|x|~ 2 99q, \x\ 2 (f(s)) ds = / (<pq,<p(s))ds. 
-oo J — oo 



In view of (I1.14p . this proves the assertion. 
Proposition 4.4. Suppose that ef G C$(R 3 \{0}). Then 



□ 



p(3) (/) 



d 3 g 



J(ci 2 -E )s 



(fi,E(0,s)a*(/)fi) ds 



4vT 2 



/ A ^ q ,x<^>- / |k| V G (k,A)/(k,A)da(k) 

JM 3 J|k|=a 2 -Eo \ i n 



\k\=q 2 -E A=12 

where dcr(k) denotes the surface measure of the sphere |k| = q 2 — Eq in 
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Proof. We start with the expression (|4.7f) for p( 3 \f) and we shall apply Lemma 14.31 to 

<p(s) = x(s)^ • (fi,E(0, *)a*(/)fi) . (4.8) 

By Lemma x(s)(p e i = e^ Hel ~ E °^ s x(p e i belongs toD(|x| 2 ) and |||x| 2 x(s)(/? e z|| < C(l + s 2 ). 
On the other hand 

(n,E(0,a)o*(/)n) = Yl I ^(k)e^ (k)s G (k,A)/(k,A)d 3 A; 

A=l,2 ^ 

= / due - *" / iu V G (k,A)/(k, A)dcr(k) (4.9) 



A=l,2 



is the Fourier transform of a function from Cg°(]R+), and hence rapidly decreasing as s —* oo. 
It follows that (|4.8p satisfies the hypotheses of Lemma 14.31 Hence Lemma 14.31 proves the 
first asserted equation because 

(^ q ,x( S )^)=e^ 2 -^ (pq.xvu). 



The second equation follows from the first one and from (|4.9p by an application of the 
Fourier inversion theorem. □ 

4.3 Proof of Theorem Q 

Lemma 4.5. 7/ / = ..,/„) u»ft e/i,...,e/ n G C °°(R 3 \{0}) and F G Cg°((-oo,0)), 

a* + (£)F(H a ) - a*(£)F(H ) = O(J) 

Proof. Choose fl6l, such that supp(/i), supp(/ n ) G {|k| < R} and then choose G G 
Co°(M) with G = 1 on supp(F) + [0, nR}. Then, by the pull through formula for a*(f) and 
by [15] . Theorem 4 (iv), 

a*(f)F(H ) = G(H )a*(f)F(H ), a\{f_)F{H a ) = G{H a )a\{f_)F(H a ). 

Using that F(Hq) — F(H a ) = O^a?), by the Helffer-Sjostrand functional calculus, that 
(a*(/) - a* + (f))F(H a ) = C(a 3 / 2 ), by the proof of Proposition and that G(H )a*(f), 
a* lr (f )F(H a ) are bounded by Lemma [A.ll and [22J Proposition 2.1, we find that 

a*(f)F(H ) - a* + {f_)F(H a ) = G(H )a*(f)F(H ) - G{H a )a\(f_)F(H a ) 
= G(H )a*(f)(F(H ) - F(H a )j + G(# )(a*(/) - a* + (f))F(H a ) 
+ (g(H ) - G(H a )y* + (l)F(H a ) = O(J) 
as a — ► 0. □ 
Recall from the introduction that is the closure of the span of all vectors of the form 

h = (h 1 ,...,h n ), where hi,whi G I^QR 3 X {1,2}), (4.10) 
and that P? is the orthogonal projection onto 
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Proof of Theorem \1.2l In the first two steps of this proof we shall establish (11.130 in the 

weak operator topology. Then we establish norm convergence to conclude the proof. 

Step 1 : Suppose H A <p = \<p, n G N and / = (/i,. ..,/„) with ef u ...,ef n E C °°(R 3 \{0}). 
Then 

lim P?(ip® a*(f)ri) = tp® a*(/)0 (4.11) 

and the analog statement holds for (p ® f2. 

Since A < there exists F £ Cq°(R) with F(X) = 1 and supp(F) C (-oo,0). Moreover 
P+F(H a ) = F(H a ) by the hypothesis of Theorem 11.21 and because S Q > for all a G R. 
Using, in addition, that 



a\(f_)F{H a ) - a*(l)F(H ) = 0(a*), 

which we know from Lemma 14.51 we conclude that 

P$(tp ® a*(/)0) = P%a*(l)F{H )(p ® Q = P%a%{f)F(H a )y ® n + 0(J) 
= a* + (l)F(H a )<p ®n + O(at) = y> ® a*(/)fi + 0(a§). 

Step 1 implies that 

lim P%$> = $ for all $ € Ranl pp (# cl ) (8) \ T . 
Step 2: to - lim P*(l c (H el ) ® l r ) = 0. 

ct^O 

Since ||Pf(l c (F e i) <g> ljr)|| < 1 for all a 6 R it suffices to show that 

lim i «(/)*«, i?(l c (fld) ® 1^) = 
q^u — 

for all p € W and all / = (/i, ...,/„) with e/i,...,e/ n G C °°(R 3 \{0}). Since o*(/)* Q G 
RanP", this follows from 

<(/)$. = a*(/)cl> o + 0( a §), 

which follows from Lemma 14.51 and Lemma IA.51 
From Step 1 and Step 2 it follows that 

u> - lim P« = l pp (H el ) \ T . (4.12) 

a— >0 

Since P? and l pp (-ff e i) (8) 1^- are orthogonal projectors, we have 

||P>|| 2 = (ip,P$ip) °^ (<P,lpp(H el ) ® l^tp) = \\l pp (H el ) ® l^f. 
Combined with (|4.12|> this proves the desired strong convergence. □ 
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5 Space-Time Analysis of the Ionization Process 

The purpose of this section is to connect our result with those of the previous papers [3 [33], 
where expressions for the zeroth and first non-trivial order of the ionization probability were 
defined. We transcribe the definitions from [33] to our model and prove their equivalence 
to the definitions in this paper. Let F R := 1{| x |>_r} ® If- 

Proposition 5.1. Let ef 1 ,...,ef N G C$(M 3 \{0}, C 3 ). Then 

lim limsup sup\\F R a*±(f )<S> a f = 0. (5.1) 

Remarks. The left hand side of Equation (|5.ip may be interpreted as the ionization 
probability to zeroth order in a [33j. Proposition 15.11 should be compared to Theorem 4.1 
in [55] . 

Proof. As in the proof of Theorem 14.11 

±oo 

a* ± (l T )^ a -a*(l T )<!> a = i J e is ( H «- E ^[W,a*(l T+s )]$ a ds, 

o 

3 

where the integral is 0{a^) in norm, uniformly in r, by Lemma 13.21 Hence it remains to 
show that 

lim lim sup sup 1 1 F R a* (/ )$ a f = 0. (5.2) 

R-+°° a \o t6 k T 

To this end, we observe that, according to Lemma lA.ll 
\\F R a*(i T )<j> a f < \\a*(f T ) 2 $ a \\\\F R $ a \\ 

N 1 

< C 2N J[ \\fiWl \\{H f + l) N t> a \\ || |x|$ a ||-. 

1=1 

This proves (|5.2p . because limsup \\(Hf + l^^H and limsup || |x|$q,|| are finite by Lemma 

a->0+ Q^0+ 

[Qland bv (lA~5l). □ 

Theorem 5.2. Let efi,...,ef N £ C\ (K 3 \{0}, C 3 ), suppose cr sc (H e i) = 0, and let r(a) = 
a -0 j or some ft g (0, |). Then 

P^(f) = lim hmsupa- 3 F R [a*_(/ r(a) )$ a - a*(/ T(a) )*«| (5.3) 



R—*oo 



a\0 



Remarks. Equation (j5.3l) is to be compared with the expression defining Q( 2 \A) in Equa- 
tion (1.9) from [33]: if we set g = a 3 / 2 and r(g) = aT^ in that equation, then Q( 2 \A) 
coincides with the right hand side of (15. 3p . 

Proof. From Proposition 13.51 we know that 



,a 3 / 2 



la* (/ , ,)$ a - a*(/ , ,)$J| < C^-T = Ca 3/2+/3 , 
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hence we may replace a*(/ , by , , )$q, for the proof of (15. 3D . From Theorem l4.ll 

we know that 

jh^limsupa^lli^fa* (/ r(a) )$ Q - <(/ t(q) )^ 



lim limsup Ili^e-^^-^^C/)!! 2 
a\0 

lim limsup HFfje-^ 1 ® 1^(/)|| 2 



where 



oo N 

*(/):= A 2p( S )^0[A(O, S ),a*(/)]fids = ^ 

" 7 1 



0j (8) rji. 



(5.4) 



(5.5) 



Explicit expressions for and ^ may be taken from the proof of Proposition 14.21 e -g-! 
rji = a*(/^)f2, but they are not needed here. From (15. 5p it follows that 



F R e~ lTH ^ 



where 



l||x||>tf } e- irffel < 



AT 

'(/>=£[ 



l{|x|>i?}e" 



!=1 



1-1; 



x||<R})e^ cl la C (^el)0Z 



<8> ?7; 



+l { || x ||> /?} e cl l pp (# cl )<^. 
By the RAGE Theorem, see [32], Satz 12.8, 

= 0, 



lim sup ||l{| x |XR}e lT cl l pp (i?ei)^| 

R—fOO T^OO 



lim ||l { | x | <i?} e ir ^l ac (^ el )^| 

T—>00 L J 



0. 



(5.6) 

(5.7) 
(5.8) 



From (jOft - lpTg]) it follows, that 

lim li 

R^oo 

by Equation (|4.6|) 



lim limsup [|Fjie- iT ^*(/)|| 2 = ||(l oc (^el) $5 VW)|| 2 = P (3) (/), 



□ 



A Uniform estimates 

In this appendix we collect estimates used in the previous sections. Most of them are well- 
known for fixed a, but this is not sufficient for us: we need estimates holding uniformly for 
a in a neighborhood of a = 0. This forces us to review some of the derivations with special 
attention to the dependences on a. 

Lemma A.l. For every N € N there is a finite constant Cm such that for all h%, hjy £ 
L 2 (M 3 x{l,2}) 

N 

\\a*(h)(H f + irf\\ < C N Y[\\hiL (A.l) 

1=1 

N 

\\a*(h 1 )---a*(h l - 1 )A(ax)a*(hi +1 )---a*(h N )(H f + l)-^\\ < C N ]l\\hj\\ u (A.2) 

3=1 
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Proof. Both, (lA.lh and (|A.2h follow from Lemma 17 in [13]. We recall that A(ax) = 
a*(G x ) + a(G x ) and we note that sup xgK 3 ||G x ||w < 00 • d 

Proposition A. 2 ([15]). .For every X < e\ there exists a constant a\ > 0, such that for all 
n G N 

sup |||x| n l(_ 00jA ](iJ Q ,)|| < 00. (A. 3) 

a<ct\ 

Proposition A. 3. There exists an a > 0, such that: 

a) For all a < a 

E a :=inia{H a ) (A.4) 

is a simple eigenvalue of H a . In the following <£ Q denotes the unique normalized 
ground state of H a whose phase is determined by ($ a ,$o) > 0. 

b) For every n E N, 

sup || |x| n $ Q || < 00. (A. 5) 

c) There exists a finite constant C, such that for all a < a and all k G M 3 \{0} 

Kk)l (l + a|k|), (A.6) 

(A.7) 
(A.8) 



||a(fc)* a || 


< 


a^C 1 - 


\Eq — E a \ 


< 




$ a - $o|| 


< 





d) For every n G N, 



sup||[iZ?- 1 ,iy a ](i/ a + i)- n+1 || < oo, (A.9) 



SUp H-ff/^a +i)" n || < OO, (A.10) 
o<o 

SUp ||JTel(-Ha + ^^ll < OO- ( A -U) 

Remark. Boundedness of [H^ 1 , H](H + i)~ n and H^(H + i)~ n has previously been 
established in [13], Lemma 5, for a class of Hamiltonians -ff that includes H a . Yet, that 
results does not imply ([A.9P and (jA.lOp . and second, its proof is much more complicated 
than the proof of (|A.9p and ([A. 10]) . because H in [13] is defined in terms of a Friedrichs' 
extension. 

Proof. That -E Q = inf a(H a ) is an eigenvalue of H a , for small a, was first shown in [2]. 
Its simplicity follows from (|A.8h . which hold for every normalized ground state vector $ a 
that satisfies the phase condition ($ a ,$o) > 0. A proof of (|A,8h may be found, e.g., in 
[15] . Proposition 19, Steps 4 and 5. A weaker form of (IA.6D is given in Lemma 20 of [15] . 
but the proof there actually shows (|A.6p . Estimate (|A.7p follows from Lemma 22 in [15] 
by choosing the infrared cutoff in this lemma larger than the UV-cutoff. Finally, (|A.5|) is a 
consequence of Proposition IA.2I and (|A.7p . 
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To prove (d) we set Rq := (Hq + z) 1 and R a := (H a +i) 1 . It is a simple exercise, using 
(TO]) and the boundedness of (H f + l^pRo, to show that \\WR \\ = C(a 3 / 2 ) as a -> 0. 
Hence we may assume that sup Q , <(5 ||Wi?o|| < 1/2 after making d smaller, if necessary. It 
follows that 

|| (#o + i)Rc\\ = 11(1 + Wi2o) _1 || < (1 - llW/Zoll)- 1 < 2 (A.12) 

for all a < a. Since H c \Rq and HfRo are bounded operators, we have thus proven statement 
(d) for n = 1, (|A,9h being trivial in this case. We now proceed by induction, assuming that 
(IA.9I) and (|A.10p hold true for all positive integers smaller or equal to a given ra > 1. To 
prove (lA.lOh for ra replaced by (n + 1) we use that 

[Hf,H a ]R^ = J2(i)^H f (W)H]- l R^ 
i=i ^ ' 

= jz(%^H } {W)R a (Hf l Rl- l -[H]-\H a ]R^ (A.13) 
1=1 ^ ' 

where sup Q<3 ||ad# (W)i2 a || < oo by (IA.12j) . by explicit formulas for ad l Hf (W) and by the 
arguments above proving that ||Wi2o|| = C(a 3/2 ). Hence \\[HJ,H a ]RZ\\ < oo follows 

from (|A.13p and from the induction hypothesis. Statement (jA.10|) with ra replaced by n + 1 
now follows from H n f +l Rl +l = (H f R a ){H]Rl) - H f R a [H],H a ]R% + \ from the induction 
hypothesis, and from (|A.9[) with n replaced by ra + 1, which we have just established. □ 

Lemma A. 4. For all l,m£ N: 



B up|i(fi>+ir$ a |i 


< 


oo, 


(A.14) 


sup|i(F / + ir<x)'<i> a || 


< 


oo, 


(A.15) 


sup p 2( 3? a 


< 


oo, 


(A.16) 


snp\\(H f + l) m H e ^ a \\ 


< 


oo, 


(A.17) 


sup|i(F / + ir(x) z P ci> a || 


< 


oo. 


(A.18) 



Proof. The statements (|A~14]) and ([ATTB]) easily follow from (|A~T0|) . (jATTTj) . and (|A~7|) . 

because Q a = (H a + i)~ n & a (E a + z) n ; note that p 2 (-ff e i + is bounded by assumption 
on V. To prove (|A. 15|) we use that 

sup \\(x) l (H f + l) m <I> Q || 2 < sup ||(x) 2/ cD a || • \\(H f + l) 2m $ a \\ 

where the right hand side is finite thanks to (1A.5|) and (|A.14[) . To prove (|A.17jl we write 

H el (H f + l) m * a = J ff el (F a + i)- 1 (^ Q + i)(F / + l) m $ Q 
= H A (H a + i)- 1 [H a ,(H f + l) m ]$ a 
+H cl (H a + i)-\H f + l) m $ a (^ Q + i). 
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The vectors [H a ,(H f + l) m ]® a and (Hf + l) m § a , and the operator H e \(H a + are 
bounded, uniformly in a < a, by (fO]> . ([00)) and (|ATLTj) . This proves (|X37|> . 

The statement (|A~T8|) follows from (|A~T5j) and <[XTl6|l after moving both p's to one 
side, and both factors (x)' to the other side of the inner product \\(Hj + l) m (x.) l p& a \\ 2 = 
{(H f + l) m {x) l p<f> a ,(H f + l) m {x) l p<!> a ). □ 



The following lemma improves upon ([AT 



Lemma A. 5. For each m £ N there is a finite constant K m , such that for all a < a 



\\(H el + i)(H f + l) m ($ a -<S> )\\<K m a2. (A.19) 

Proof. Let A := (E + ei)/2. Thanks to (lA~7l) in Proposition El we may assume that 
sup Q , <5 E a < A by making a smaller, if necessary. Pick g £ C^°(W) with suppg C (— oo, A) 
and with g{E a ) = 1 for all a < a. On the one hand, 



(H cl + i){H f + l) m g(H )(* a - $ )|| 
< \\(H cl + i)(H f + irg(H )\\\\^ a -^ \\ = O(a 



3/2 ^ 



by (1A.8|) . On the other hand, (1 - g(H )){<& a - <J> ) = (g{H a ) - g(H ))$ a by construction 
of 5. Hence it remains to prove that 

\\(H el + i)(H f + l) m (g(H a ) - g(H ))$ a \\ = 0(a 3 / 2 ). (A.20) 

To do so, we use the Helffer-Sjostrand functional calculus with a compactly supported 
almost analytic extension g of g that satisfies an estimate |dz£f(z)| < C|y| 2 . Here and 
henceforth z = x + iy with x,y £M. It follows that 

(H el + i){H f + l) m (g(H a ) - g(H ))^> a 
= -- f (H el + i)(H -z)- l (H f + l) m W(H a -zy l <S> a %dxdy (A.21) 

where 

(H f + l) m W = Y J [ r ] l )ad l Hf (W)(H f + l) m - 1 =: a 3 / 2 W(m)(H f + l) m . (A.22) 
1=0 ^ 1 J 

From the equations [Hf, a*(G x )} = a*(uG x ) and [Hf,a(G x )] = —a(ujG x ) it is clear that 
the operator W(m), defined by (|A.22p . is i^o-bounded. Hence we can estimate the norm of 
(|A.2ip from above by 



a 3 / 2 



J/2 r- 

— \\(H cl + i)(H + i)- 1 !! J 



Oz 



H + i 



H -z 



\z-E n 



dxdy 



x \\W(m)(H +i)- l \\\\{H + i)(H f + l) m <f> a \\. (A.23) 



The integral is finite by construction of g, because \z — E^ 1 < |y| _1 , and because \\(Hq + 
i)[Hq — z)~ 1 \\ < 1 + (1 + |x|)/|y| by the spectral theorem. The last factor in (jA.23|) is 
bounded uniformly in a < a by (|A. 14|) and (|A.17|) from Lemma IA.41 This establishes 
(|A.20p and thus concludes the proof of the lemma. □ 
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Lemma A. 6. Suppose that V satisfies the hypotheses in Section^ Then 

(i) xip e i G D(H cl ) and (i? el - E )xip e i = -2Vip el . 

(ii) e~ lHclt ~x.<p e i G -D(|x| 2 ) and there exists a constant C such that for all t G M ; 

||| x |2 e -«rei* x ^|| < C(l+t 2 ). 

Proof, (i) For all 7 G C^°(M 3 ) we have x# e i7 = H e \x.^ + 2V7 and hence 

= ('y,E xip e i - 2V(f e l) . 

Since Cq°(IR 3 ) is a core of H e \, we conclude that x<p e ; G D(H C \) and that 

i? e lX</? eZ = E Xif e i - 2V</?. 

(ii) Let V> := £j<£eZ for some i G {1, 2, 3}. We shall only need that ift G L>(|x| 2 ) n D(-A) 
which follows from (i). By the fundamental theorem of calculus, in a weak sense 

t 



t 

= |x|fy + 2^ e isHol ( x • P + P • x) e - isHol ^s 


= |x| 2 ^ + 2t(x • p + p • x)il> + 2 J ds j dre irH ^(4p 2 - x • VV)e~ irHcl ilj. (A.24) 



Here ip G L>(|x| 2 ) n D(-A) C L>(x • p + p • x) and e~ irH ^ G D(fT d ) = D(-A) because 
■0 G D(H e \) by part (i). Therefore assertion (ii) follows from (|A.24p and from the hypotheses 
on V. □ 
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